GT strengths and electron-capture rates for p/-shell nuclei of relevance for late stellar 

evolution. 



A.L. ColcE and T.S. Anderson 
Physics Department, Kalamazoo College, Kalamazoo, MI 49006, USA 

R.G.T. ZegersQ Sam M. Austin, B.A. Brown, and L. Valdez 
National Superconducting Cyclotron Laboratory, Michigan State University, East Lansing, MI 48824, USA 
Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA and 
Joint Institute for Nuclear Astrophysics, Michigan State University, East Lansing, MI 48824, USA 

S. GuptfH 

Theoretical Division, Los Alamos National Laboratory, Los Alamos, NM 87545 

G.W. Hitt and 0. Fawwaz 
Department of Applied Mathematics and Sciences, Khalifa University of Science, 
Technology, and Research, P.O. Box 127788 Abu Dhabi, UAE 
(Dated: March 2, 2013) 

Background: Electron-capture reaction rates on medium-heavy nuclei are an important ingredient for modeling 
the late evolution of stars that become core-collapse or thermonuclear supernovae. The estimation of these rates 
requires the knowledge of Gamow- Teller strength distributions in the /3 + direction. Astrophysical models rely 
on electron-capture rate tables largely based on theoretical models, which must be validated and tested against 
experimental results. 

Purpose: This paper presents a systematic evaluation of the ability of theoretical models to reproduce exper- 
imental Gamow- Teller transition strength distributions measured via (n,p)-type charge-exchange reactions at 
intermediate beam energies. The focus is on transitions from stable nuclei in the pf shell (45 < A < 64). In 
addition, the impact of deviations between experimental and theoretical Gamow- Teller strength distributions on 
derived stellar electron-capture rates is investigated. 

Method: Data on Gamow- Teller transitions from 13 nuclei in the pf shell measured via charge-exchange reactions 
and supplemented with results from /3-decay experiments where available, were compiled and compared with 
strength distributions calculated in shell-models (using the GXPFla and KB3G effective interactions) and quasi- 
particle random-phase approximation (QRPA) using ground-state deformation parameters and masses from the 
finite-range droplet model. Electron-capture rates at relevant stellar temperatures and densities were derived for 
all distributions and compared. 

Results: With few exceptions, shell-model calculations in the pf model space with the KB3G and GXPFla 
interactions qualitatively reproduce experimental Gamow- Teller strength distributions of 13 stable isotopes with 
45 < A < 64. Results from QRPA calculations exhibit much larger deviations from the data and overestimate the 
total experimental Gamow- Teller strengths. For stellar densities in excess of 10 7 g/cm 3 , ground-state electron- 
capture rates derived from the shell-model calculations using the KB3G (GXPFla) interaction deviate on average 
less than 47% (31%) from those derived from experimental data for which the location of daughter states at 
low excitation energies are well established. For electron-capture rates derived from Gamow- Teller strengths 
calculated in QRPA, the deviations are much larger, especially at low stellar densities. 

Conclusions: Based on the limited set of test cases available for nuclei in the pf shell, shell-models using the 
GXPFla and KB3G interactions can be used to estimate electron-capture rates for astrophysical purposes with 
relatively good accuracy. Measures of the uncertainties in these rates can serve as input for sensitivity studies 
in stellar evolution models. Ground-state electron-capture rates based on the QRPA formalism discussed in the 
paper exhibit much larger deviations than those based on the shell-model calculations and should be used with 
caution, especially at low stellar densities. 

PACS numbers: 21.60.Cs, 21.60.Jz, 23.40.Hc, 25.40.Kv, 26.30.Jk 
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Supernovae play a critical role in the universe. They 
leave behind neutron stars and black holes and their 
Shockwaves drive mixing of interstellar material and stim- 
ulate galactic chemical evolution. Supernovae are also 
major sources of nucleosynthesis. However, for the two 
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main types of supernovae, core-collapse (type II) and 
thermonuclear (type la) supernovae, the driving mech- 
anisms are still not well understood. A variety of nu- 
clear physics processes play important roles in both types 
and improving the nuclear physics input for astrophysical 
simulations of supernovae has been a strong motivation 
for a wide variety of nuclear experiments and theoretical 
calculations. 

Electron-captures (EC) on medium-heavy nuclei play 
an important role in both types of supernovae [l[ . During 
the pre-explosion evolution of core-collapse supernovae, 
when the Fermi energy of the degenerate electron gas be- 
comes sufficiently high to overcome Q-value limitations 
that restrict EC under terrestrial conditions, the nuclear 
matter in the stellar core is neutronized and the electron 
abundance (Y e ) is reduced by EC reactions [2|]. Conse- 
quently, the electron pressure is reduced, and the energy 
and entropy drop due to the emission of neutrinos in the 
EC reactions [H, H, H| . In the early pre-explosion evolu- 
tion, EC on nuclei in the pf shell are key. Just prior to 
the collapse, EC on nuclei in the sdg shell also plays an 
important role. 

In type la supernovae, which are thought to occur in 
binary systems in which a white dwarf accretes material 
from a companion star, the thermonuclear flame leaves 
behind an equilibrium distribution of iron-region (jpf- 
shell) nuclei. Electron-captures on these nuclei reduce 
the pressure and retard the expansion of the star. EC 
reactions also reduce the amount of 56 Ni produced, mak- 
ing the explosion less bright (see e.g. [1], and references 
therein). Type la supernovae are thought to produce 
about half of the iron-group nuclei in the solar system. 
Hence, the nature of their progenitors and the explosion 
models can be constrained by the condition that these su- 
pernovae should not generate amounts of these elements 
in excess of a factor of two compared to solar system 
abundances @. However, to reliably do so, it is key to 
use accurate weak-reaction rates in the simulations [5j . 

To establish an accurate database of weak reaction 
rates of importance for astrophysical simulations, re- 
liable Gamow- Teller (GT) transition-strength distribu- 
tions must be available for a large number of nuclei. Ad- 
ditionally, EC in stellar environments occurs at relatively 
high stellar temperatures and transitions from thermally 
populated excited states can contribute significantly and 
must also be taken into account. It will, therefore, be im- 
possible to measure a large fraction of the relevant GT 
strength distributions and one has to rely on theoretical 
estimates. 

In the past, astrophysical simulations often used EC 
rates calculated by Fuller, Fowler and Newman (FFN) 0- 
[Toj . These rates were based on GT strengths calculated 
in an independent particle model (IPM), supplemented 
with experimental strengths from /3-decay experiments. 
Based on results from charge-exchange experiments (see 
e.g. Ref. 11]), it has become evident that GT tran- 
sition strengths are strongly fragmented and quenched 
compared to the IPM calculations, largely because of the 



residual interactions between the nuclear constituents. 
Modern calculations take into account such residual in- 
teractions and updated libraries for weak-reaction rates 
based on shell-model calculations have been produced for 
nuclei in the sd [l2[ and pf shells [H| . For nuclei beyond 
the pf shell, schemes beyond the shell-model are required 
(see e.g. [jj, and references therein, for an overview). 

In this work, we focus on the accuracy of estimates for 
EC rates on a variety of p/-shell nuclei (45 < A < 64) 
based on theoretical GT strength distributions. One can 
measure GT strength distributions directly through j3- 
decay experiments, but the excitation-energy range that 
can be covered is severely limited by the available Q- 
value window. Therefore, charge-exchange reactions at 
intermediate energies (~ 100 MeV/u and above) have 
become the preferred tool to probe GT strength distribu- 
tions over the full range of excitation energies of interest. 
The charge- exchange experiments rely on the proportion- 
ality between GT transition strengths and the measured 
differential cross sections extrapolated to values at van- 
ishing linear momentum transfer ([3]; see also Section 
HI)) . For testing theoretical calculations of GT strength 
distributions in p/-shell nuclei of relevance for EC rates 
in late stellar evolution, measurements have been per- 
formed using the (n,p), (d, 2 He) and (i, 3 He) reactions. 

Although comparisons of GT strengths and EC rates 
derived from data and theory have been made for spe- 
cific nuclei or for subgroups of nuclei with specific the- 
oretical models in the past, a comprehensive evaluation 
for nuclei in the pf shell is lacking. Such an evalua- 
tion is important to gain better insight in the deficien- 
cies of theoretical models, to obtain a measure of the 
uncertainties in the EC rates used in astrophysical sim- 
ulations, and to provide guidance to future experimental 
and theoretical work. The goals of the current work are 
two-fold: (i) to perform a consistent comparison of the 
experimentally extracted GT strength distribution from 
charge-exchange reactions on stable p/-shell nuclei with 
results from various theoretical models, and (ii) to study 
the impact of differences between the theoretical and ex- 
perimental strength distributions on the estimates for EC 
rates used in astrophysical modeling. Besides the data 
from the above-mentioned AT Z = +1 charge-exchange 
probes on stable p/-shell nuclei, we also include data from 
AT 2 = — 1 (p,n) charge-exchange experiments. Relying 
on the fact that isospin-symmetry breaking effects are 
small, they can also be used to extract GT strength dis- 
tributions of relevance for estimating EC rates. Where 
available, we also included /3-decay data. In total, 25 dis- 
tinct data sets were used, for 13 separate nuclei. These 
nuclei and the data set type(s) are listed in Table H 

Three sets of theoretical calculations have been tested 
against experiment in the present work. Two of those sets 
are based on shell-model calculations in the pf model 
space. These shell-model calculations were performed 
with the code NuShellX@MSU [11, OJ]. One set of cal- 
culations employed the GXPFla interaction [l7l - fl9l |. the 
other the KB3G interaction 1201. The latter interaction 
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is an updated version of the KBF interaction [2l| , which 
was used to generate the weak reaction rate library of 
Ref. [l3[ and whose parameters were primarily deduced 
from experimental data in the lower pf shell. Parameters 
of the GXPFla interaction have been fitted to reproduce 
the experimental excitation energies and masses for many 
p/-shell nuclei. 

The third set of theoretical GT strength distributions 
was based on the Quasi-particle Random Phase Approx- 
imation (QRPA) formalism of Ref. 22], using ground- 
state deformation parameters and masses from the finite- 
range droplet model of Ref. [231 ■ strengths produced 
in this model have been used, for example, to estimate 
heating in the accreted neutron star ocean [24|. In the 
following, we refer to this set of calculations as "QRPA" . 
Although not the focus of the current work, usage of the 
this QRPA formalism for calculating GT strength distri- 
butions has the advantage that it can be applied to nuclei 
beyond the p/-shell. 

This paper is organized as follows. In Section HH we 
briefly review the extraction of GT strengths from the 
data and the manner in which the experimental and the- 
oretical GT strength distributions were entered in the 
EC-rate calculation. The EC rate calculations are also 
briefly described, focusing on input parameters related to 
the GT strength distributions. In Section Unl a case- by- 
case and self-contained comparison of experimental and 
theoretical GT strengths is provided for each of the nuclei 
listed in Table HI In addition, derived EC rates are also 
calculated for different astrophysical conditions. In Sec- 
tion llVl the results from the comparisons between theory 
and experiments for individual cases are combined in an 
attempt to gain a better insight into the overall capabil- 
ity of the theoretical models to accurately estimate EC 
rates of astrophysical interest. 



II. GT STRENGTHS AND 
ELECTRON-CAPTURE RATES 

If data are available for the decay half-life, GT transi- 
tion strengths (B(GT)) can be calculated directly from 
the comparative half-life ft using: 



cross section, extrapolated to zero linear momentum 
transfer (q = 0) |14| : 



(^-) 2 B(GT) 
9v 



K/g 2 v 

ft : 



(1) 



where ^ = -1.2694 ± 0.0028 [2j and Kjg 2 v = 6143 ±2 
s j2(3 . Here, B{GT) is defined such that it equals 3 for 
the decay of the free neutron. Of the cases studied here 
(see Table IJ), (3— decay data were only available for four 
nuclei and only provided information about the transition 
between the ground states of the mother and daughter 
nuclei. 

As mentioned above, the extraction of GT strengths 
from intermediate-energy (E > 100 MeV/u) charge- 
exchange reactions relies on the proportionality between 
the transition strength and the measured differential 



= &B(GT), 



(2) 



GT 



where a is the unit cross section. The unit cross sec- 
tion can be derived from direct comparison of the charge- 
exchange cross section to the B(GT) extracted from /3- 
decay data, if the latter are available. Otherwise, one 
relies on the unit cross section extracted from a nucleus 
of similar mass number, or an empirically established 
correlation between mass number and unit cross section 

[II SMI. 

When the excitation-energy resolution in the charge- 
exchange experiments is high compared to the level den- 
sity, GT strengths can be extracted for each transition to 
a specific final state. In general, however, this in not the 
case for the experimental data considered here, or is only 
true at low excitation energies. Therefore, the GT contri- 
butions to the excitation-energy spectra were extracted 
using a Multipole-Decomposition Analysis (MDA); the 
AL — 0, AS" = 1 GT component of the total response 
in a given excitation-energy bin or peak was determined 
by fitting the angular distribution to a linear sum of the- 
oretically calculated angular distributions, each associ- 
ated with different angular momentum transfer. Conse- 
quently, the exact location of the GT strength is not nec- 
essarily known accurately due to the limited experimen- 
tal excitation-energy resolution and the analysis of the 
data into excitation-energy bins of finite width. When 
preparing the GT transition strengths for the calcula- 
tion of EC rates, we used the same energy bin widths 
as used in the original analyses of the data if transitions 
to specific final states could not be isolated and placed 
all strengths in the center of the bins. The GT strength 
distributions are also represented in this manner in the 
figures. Exceptions were made for the cases where, due 
to the limited energy resolution, GT strength appeared 
below the ground state of the final nucleus. Details are 
provided on a case-by-case basis in Section IIIII 

Besides the uncertainty in the unit cross section (ap- 
proximately 10% |l4 l28l430l |). typical uncertainties in the 
extraction of the GT strength for a transition to a par- 
ticular final state are 10 — 20%, with error bars increas- 
ing for weaker transitions strengths [H], HH |32{ . The 
dominant source of these uncertainties is related to the 
interference between AL = and AL = 2 amplitudes 
which both can contribute to the A J = 1 GT transitions. 
The interference between these amplitudes is caused by 
the tensor component of the effective nucleon-nucleon 
interaction that mediates the charge-exchange reaction 
[Til [33| . Since the interference can be constructive as 
well as destructive, depending on the exact nature of 
the wave functions of the initial and final states, these 
uncertainties in the extraction of the strength for spe- 
cific states tend to cancel when summing over several 
states [3l], HH ■ The uncertainties in the magnitudes of 
the extracted GT strengths are not explicitly referred 
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TABLE I. Overview of the GT transitions studied in this work. For each of the cases, the data types available and the types 
of theoretical calculations used are indicated. In addition, the ground-state to ground-state Q-value for electron-capture is 
provided. The last column gives the references of the papers from which the data sets used in the comparisons to theory were 
drawn. 
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X 
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[61-64] 


64 Zn(0+) 


64 Cu(l+) 


X 
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X 




X 


X 
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-1.09 


[32, 65, 66] 



a Using T> transitions and applying isospin symmetry (see text). 

b Shell-model calculations were performed in truncated model space (see text). 

c Qec is calculated using nuclear masses. 



to in the remainder of this paper. Since EC rates are 
proportional to GT strengths (see below) , the uncertain- 
ties in the GT strengths extracted from the data limit a 
meaningful comparison between EC rates based on the- 
oretical and experimental GT strength distribution to a 
10-20% level of discrepancy. As discussed below, the EC 
rate is much more sensitive to the excitation energies of 
the states in the daughter nuclei to which EC takes place. 

In the shell-model calculations of the GT strengths of 
relevance for this work, transitions to about 100 final 
states of a given spin-parity were generated, which was 
sufficient to cover the excitation energy ranges for which 
experimental data were available and to ensure that the 
great majority of the total strength (always more than 
95%) was included in the theoretical results. Except for 
two cases (detailed below) for which the computational 
time would have been excessive, the shell-model calcula- 
tions were performed in the full pf shell, i.e. assuming a 
40 Ca inert core. The transition to the ground state of the 
daughter nucleus is not always of GT nature. Therefore, 
the level scheme of each daughter nucleus at low excita- 
tion energies was also generated in the shell-model calcu- 
lations (separately for each interaction) to ensure a con- 
sistent placement of the first level associated with a GT 
transition from the parent relative to the ground state 
calculated within the model. No attempt was made to 
shift the theoretical strength distribution to better match 
the data since this cannot be done for the great majority 
of transitions relevant to astrophysical applications for 
which no data are available and one thus has to rely on 
theory only. The calculated strengths in the shell-models 
were scaled by a quenching factor which was estimated 
at (0.74) 2 for nuclei in the pf shell [34j. The quench- 



ing factor accounts for degrees of freedom not included 
in the shell-model calculations. A large fraction of that 
quenched strength is known to be located at higher exci- 
tation energies, due to the mixing between lp — lh and 
2p — 2h configurations (see e.g . Ref. [35[ and references 
therein). Following Ref. [22|], a quenching factor was 
not applied to the GT strengths calculated in the QRPA 
framework; these strength distributions were not modi- 
fied in any fashion prior to the comparison with the data 
or the calculation of the EC rates. 

The EC rate calculations were performed following the 
formalism described by Fuller, Fowler and Newman 0- 
[loj . which was implemented in a code previously used 
in Refs. 0, We briefly review the formalism for 
calculating EC rates of Refs. 0413; insofar as they are 
helpful for the further discussion in this paper. 

The EC rate (A) is calculated with: 



A = ln2^ 



f 3 (T,P,U F ) 
(ft)j ' 



(3) 



where the sum runs over all daughter states that can be 
populated through a GT transition in the EC reaction. 
Only transitions from the ground state of the parent nu- 
cleus are considered in the calculations presented in this 
paper. The phase space integral fj depends on the tem- 
perature T, density p and electrochemical potential Uf- 
It is calculated by: 



W- 2 



Qec — Ej 
m e c 2 



+ w) G(Z,w)fedw, (4) 



where Qec is the ground-state to ground-state EC Q- 
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value (calculated with nuclear masses), Ej is excitation 
energy of the j th daughter state, m e c 2 is the electron rest 
mass, w = E e /m e c 2 the total electron energy in units of 
m e c 2 , and G(Z, w) is a factor that accounts for the effects 
of the Coulomb barrier. The lower limit wi corresponds 
to the threshold value of w for which EC becomes ener- 
getically possible. Realizing that the blocking in the final 
neutrino phase space plays a negligible role, and neglect- 
ing any effects related to the presence of bound electrons 
and ions, f e represents the Fermi-Dirac distribution for 
the energies of the electrons in the stellar plasma: 



fe 



exp 



U-U f 



1 



(5) 



U is the kinetic energy of the electron, and Up/ksT is 
referred to as the degeneracy parameter, and ks is the 
Boltzmann constant. At T = 0, 



U F (T = 0) = 0.511 (l.018(p 6 ^e) 3 +1) -1 



(6) 



where pe is the density divided by 10 6 g/cm 3 and Y e is 
the electron fraction. If the Fermi energy Ep — Up+m e c 2 
exceeds Qr = E% — Qpc, where E\ is the excitation en- 
ergy of the first state that can be captured into in the 
daughter nucleus, EC can occur. Eq. © shows that with 
increasing stellar density, capture to daughter states at 
higher excitation energies becomes increasingly impor- 
tant. With increasing temperature (decreasing degener- 
acy parameters), the Fermi surface smears out, and EC 
can occur even when e^(T=0) < Qr. If that is the case, 
the EC rate is particularly sensitive to the excitation en- 
ergy of the daughter states, and small differences between 
experimentally and theoretically deduced excitation en- 
ergies can lead to very large differences in derived EC 
rates. Therefore, the comparison between EC rates de- 
rived from data and theory is difficult if the experimental 
excitation-energy resolution is poor. Similarly, the abil- 
ity of a particular theoretical model to accurately predict 
the excitation energies of the daughter states is critical 
for the estimation of accurate EC rates. 

For the current work, EC rates were calculated over 
a wide stellar density-temperature grid, where T var- 
ied between 10 7 K to 10 11 K and pY e varied between 
10 and 10 14 g/cm 3 . However, in the following, we fo- 
cus on two particular densities: pY e = 10 7 g/cm 3 (case 
I) and 10 9 g/cm 3 (case II), which roughly delineate the 
density range of interest for EC rates in supernovae. Case 
I (at T w 3x 10 9 K) is representative for the conditions of 
the core during the silicon-burning phase of a pre-core- 
collape progenitor star |3jj. Case II (at T ~ 10xl0 9 K) 
is representative for the conditions present just prior to 
the collapse of the core [1, [H| . Case II is is also repre- 
sentative for the high-density burning regions in which 
EC occurs during the thermonuclear runaway in Type la 
supernovae [!, Q . 

At pY e = 10 7 (10 9 ) g/cm 3 , e F (T=0) » 1.2(5.2) MeV. 
For many of the nuclei discussed in the paper, Qr > ep 
at the lower of these densities and the EC rates are very 



sensitive to captures to states at low excitation energies in 
the daughter nuclei and the temperature. At the higher 
density, capture into states at higher excitation energies 
can, depending on the value of Qpc, play a significant 
role. 



III. COMPARISON OF EXPERIMENTAL AND 
THEORETICAL GT STRENGTH 
DISTRIBUTIONS 

In this section, we compare the experimental and the- 
oretical GT strength distributions, and the derived EC 
rates, for each of the cases listed in Table |U 



>Sc- 



3 Ca 



Fig. [TJ shows the experimental and theoretical GT 
strength distribution for transitions from the 45 Sc(^ ) 
ground state. States in 45 Ca with J n = (|, |, |)~ 
are populated. The B(GT) for the transition to the 

45 Ca(| ) ground state is known from /3-decay (4.9 x 
10~ 3 ) [33]. The GT strength distribution has been mea- 
sured in a 45 Sc(n,p) experiment [40| . Because of the 
limited resolution of that experiment, results are avail- 
able only in excitation-energy bins of 1 MeV (0.5 MeV 
for the energy bin between i?^( 45 Ca)=0-0.5 MeV). In 
Ref. 40] it was noted that, due to the possible con- 
tamination from the (n,p) reaction on hydrogen in the 
target, the extracted strength up to E x ( ib G&)=2 MeV 
is an upper limit. For the present study, the data avail- 
able from /3-decay and 45 Sc(n,p) were combined. GT 
strength in excess of the j3— decay strength reported be- 
tween E x ( 45 Ca) =0-0.5 MeV in the 45 Sc(n,p) experiment 
was placed at £ x ( 45 Ca)=0.25 MeV. For E X ( 45 C&)> 0.5 
MeV, the B(GT) values from the 45 Sc(n,p) experiment 
were used, placed at excitation energies in the center of 
each bin. The combined plot is shown in Fig. [TJa). 

Fig. \V[b) shows the results from the QRPA calculation, 
divided by a factor of 3. Such rescaling of the QRPA 
strength distributions (for visualization purposes only) 
is performed throughout the paper. Figs. QJc) and (d) 
show the results from the shell-model calculations in the 
full pf shell-model space using the KB3G and GXPFla 
interactions, respectively. Fig. [IJe) shows the running 
sum of B(GT) as a function of excitation energy for the 
data and each of the theoretical calculations. 

Under the assumption that the broad distribution seen 
around 7 MeV in the (n,p) data is from a transition to 
a single state in 45 Ca (but broadened due to the limited 
excitation energy resolution), all three theoretical calcu- 
lations predict the location of that state about 1 MeV 
too low. About 0.2 units of GT strength are found be- 
low 2 MeV in the experiment. The shell-model calcula- 
tions using the GXPF1A (KB3G) interaction predict a 
summed B(GT) of 0.014 (0.033) below E x { 45 Ca)=5 MeV 
and a B(GT) of 0.00023 (0.00019) for the ground-state 
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E x rCa) (MeV) 



FIG. 1. GT transitions from 45 Sc to 45 Ca: (a) from (n,p) 
and p— decay data, (b) from calculations in QRPA (divided 
by a factor of 3), (c) from shell-model calculation using the 
KB3G interaction and (d) from shell-model calculations using 
the GXPFIa interaction. In (e), running sums of B(GT) as a 
function of excitation energy are plotted. Here, and all similar 
graphs, the shell-model strengths have been scaled by (0.74) 2 
(see text). 



to ground-state transition, which is about a factor of 10 
smaller than the experimental value and hard to distin- 
guish in Figs, [ljc) and (d). Above 5 MeV, much better 
consistency is found, although the summed strength up 
to 10 MeV is better reproduced in the shell- model cal- 
culation using the GXPFIa interaction than the calcu- 
lation using the KB3G interaction. In the QRPA cal- 
culations, the ground-state to ground-state transition is 
about thirty times stronger than deduced from the /3- 
decay measurement and the total GT strength is overes- 
timated by a factor of 3. 

Fig. (Ha) and (b) show the EC rates on the ground 
state of 45 Sc as function of stellar temperature, based on 
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"— data 45 Sc(n,p)+ 45 Ca(P") 


GXPFIa 


; KB3G 


QRPA 



Temperature(10 K) 

FIG. 2. Electron-capture rates for 45 Sc(e~, z/ e ) 45 Ca based on 
experimental and theoretical GT strength distribution from 
Fig. [TJ plotted as a function of temperature at pY e = 10 7 
g/cm 3 (a) and pY e =10 9 g/cm 3 (b). 



the three theoretical and the experimental GT strength 
distributions shown in Fig. [TJ at values of pY e of 10 7 
g/cm 3 and 10 9 g/cm 3 , respectively. Because Qec=- 
0.7677 MeV (See Table[T|, states in 45 Ca at low excitation 
energies reside below the Fermi surface (sp(T = 0) = 1.2 
MeV) even at the lower density, and the derived EC rates 
do not change drastically as a function of temperature. 
However, the excess of strength seen at low excitation 
energies in the experimental distribution gives rise to a 
strongly enhanced (x20 — 100) EC rate in comparison 
to those derived from the shell-model calculations. Con- 
versely,the strong transition to the ground state of 45 Ca 
in the case of the QRPA calculation results in a much 
higher rate. At the higher density (Fig. [^b)) the transi- 
tions to the excited states start to play a role, including 
the strong transition to the state at E x ( A5 Sc)^ 6 MeV. 
Nevertheless, since the phase-space is small for a tran- 
sition at such a high excitation energy the discrepancies 
between the EC rates seen at the lower density remain at 
the higher density. Compared to the other nuclei studied, 
the case of 45 Sc exhibits amongst the largest discrepan- 
cies between the data and shell-model calculations. Early 
work [4l| indicated that mixing between pf and sd-shell 
configurations can occur for states at relatively low exci- 
tation energies in the lightest pj-shell nuclei, which could 
affect the GT strength distributions. To confirm and 
evaluate the impact of possible configuration mixing, fur- 
ther experimental data would be helpful. In addition, it 
is important to resolve the possible contamination of the 
experimental GT strength distribution at low-excitation 
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energies from reactions on hydrogen in the 45 Sc target. 



Fig. [3] shows the experimental and theoretical GT 
strength distributions for transitions from the 48 Ti(0 + ) 
ground state to 1 + final states in 48 Sc. GT strengths 
from two sets of 48 Ti(n,p) data are available. In the ear- 
lier experiment !42j performed at E(n) — 200 MeV, the 
excitation-energy resolution was 1.2 MeV. In the later ex- 
periment [43j], performed at E(n) — 300 MeV, a similar 
resolution was achieved, but with better statistical accu- 
racy. In a comparison of the two data sets, we found that 
below E x ( 48 Sc)=6 MeV, the GT strength distributions 
were consistent, but that at higher excitation energies, 
slightly less strength was reported in Ref. [43j]. Since 
a higher statistical accuracy allows for a more accurate 
MDA, we used the data (and adopted the same energy 
binning) from Ref. [43] in our analysis, which is shown 
in Fig. EJa). 

Fig. EJb) shows the GT strength distribution extracted 
from a 48 Ti(d, 2 He) experiment at £(d)=183 MeV Q. 
The excitation-energy resolution was 120 keV, which 
made it possible to extract the GT strengths on a state- 
by-state basis up to E x ( 48 Sc)—5 MeV. GT strengths were 
not extracted beyond that energy. Since Qec — —4.505 
MeV, this limit didn't affect the EC rate calculations 
significantly, as the levels populated beyond E x ( 48 )Sc=5 
MeV are located above the Fermi level {ef{T = 0) = 5.2 
MeV), even at pY e — 10 9 g/cm 3 ) and a significant 
amount of strength is present below the Fermi level. 

Fig. Etc) shows the results from the QRPA calcula- 
tion, divided by a factor of 6. Figs. Hfd) and (e) show 
the results from the shell-model calculations in the full 
pf shell-model space using the KB3G and GXPFla in- 
teractions, respectively. Fig. EJf) shows the running sum 
of B(GT) as a function of excitation energy. 

Both shell-model calculations do about equally well in 
describing the high-resolution (<i, 2 He) data, but the cal- 
culation with the GXPFla interaction provides a better 
estimate for the total strength found up to 10 MeV in the 
(n,p) data. This is mostly due to the fact that slightly 
more strength is predicted at low excitation energies in 
the calculation with the GXPFla interaction. Neither 
the GXPFla or KB3G interaction describe the long tail 
seen in the strength distribution extracted from the (n,p) 
data. In the QRPA calculations, the strength distribu- 
tion is dominated by a single transition to a state at 5 
MeV, in clear contradiction to the data. A weaker tran- 
sition at about 1 MeV is also predicted, but is not seen 
in the data. 

The EC rates derived from the experimental and the- 
oretical strength distributions is displayed in Figs. |U 
The rates calculated from the strength distributions gen- 
erated with the GXPFla and KB3G interactions are 
similar and match well with those based on the high- 
resolution 48 Ti(d, 2 He) data. The rates based on GT 



strengths extracted from the 48 Ti(n,p) data are much 
higher, largely due to the poor energy resolution. In 
combination with the chosen binning of the data, this 
causes some of the GT strength to be placed at artificially 
low excitation energies. The rates calculated based on 
the strength distribution in QRPA are also higher than 
the rates based on the shell-model calculations and the 
48 Ti(d, 2 He) data due to the relatively strong transition 
to the above-mentioned state at about 1 MeV. There- 
fore, the correspondence between the rates based on the 
QRPA calculations and the 48 Ti(n,p) data is coinciden- 
tal. 



C. 50 V— > 50 Ti 

Information about GT transition strength from the 
50 V(6+) ground state to (5, 6, 7)+ states in 50 Ti are only 
available from a 50 V(<i, 2 He) experiment performed at 171 
MeV [45[. The excitation-energy resolution was 160 keV. 
The extracted GT strength distribution is shown in Fig. 
[5ja). A known 6 + state at 3.198 MeV [46j . was obscured 
by a strong peak from the H((i, 2 tle) reaction at forward 
scattering angles, but the authors could determine that 
its B(GT) was "insignificant" from the data at larger 
scattering angles. The location of this state is indicated 
in Fig. [5ja). Since the high level density made it im- 
possible to assign GT strengths on a state-by-state basis, 
the strengths shown are in 50-keV wide bins, following 
Fig. 4(a) from Ref. j4j|. 

Figs. Ob,c,d) shows the theoretical calculations for the 
GT strength distribution from the QRPA (divided by 3) 
and shell-model calculations in full pf model space using 
the KB3G and GXPFla interactions, respectively. Fig. 
03e) shows the running sum of B(GT) as a function of 
excitation energy. 

The shell-model calculations with the GXPFla and 
KB3G interactions are very similar and both correspond 
quite well to the experimental GT strength distribution. 
In addition, both predict a weak transition to a 6 + state 
situated about 1 MeV below transitions to other states, 
consistent with the experiment. In both shell-model cal- 
culations, somewhat less fragmentation of strength com- 
pared to the experimental results is predicted up to an ex- 
citation energy of 8 MeV. A similar, but enhanced effect 
is seen over the excitation energy range up to 11 MeV for 
the calculations performed in QRPA. Compared to the 
shell-model calculations and the data, the QRPA calcu- 
lations also predict a relatively large amount of strength 
at excitation energies above 9 MeV. 

In Figs. EJa) and (b), EC rates based on these 
strengths are displayed. Because the transition to the 
6 + state at 3.198 MeV could not be included in the 
rate based on the 50 V(<i, 2 He) experiment, the corre- 
sponding rate estimate is lower. If an EC transition to 
this 6 + state, with a strength similar to that calculated 
with the KB3G and GXPFla interactions, was artifi- 
cially included when generating the EC rates based on 
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FIG. 3. GT transitions from 48 Ti to 48 Sc: (a) from (n,p), (b) 
from (<i, 2 He) data, (c) from calculations in QRPA (divided 
by a factor of 6), (d) from shell-model calculations using the 
KB3G interaction and (e) from shell-model calculation using 
the GXPFIa interaction. In (f), running sums of B(GT) as a 
function of excitation energy are plotted. 



the 50 V((i, 2 He) experiment, the rate at pY e =10 7 g/cm 3 
matches about equally well to either of these two theo- 
retical rates. At pY e =10 7 g/cm 3 , the calculation that 
employs the GXPFIa interaction does slightly better, in- 
dependent of whether a small contribution from transi- 
tions to the first 6 + state was artificially included for the 
rate calculation based on the 50 V(e?, 2 He) data. 
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FIG. 4. Electron-capture rates for 48 Ti(e~, ^ e ) 48 Sc, based on 
experimental and theoretical GT strength distributions shown 
in Fig. [3l plotted as a function of temperature at two den- 
sities, pY e =10 7 g/cm 3 (a) and pY e =10 9 g/cm 3 (b). Note 
that the rates based on the shell-model calculations with the 
KB3G and GXPFIa interactions overlap. 



The value of Qec for the transition from the 50 V 
ground state to the 50 Ti ground state is positive, unlike 
the other cases studied here (see TablelJ). Combined with 
the weakness of the transition to the first 6 + state, the 
transitions to the higher-lying excited states thus play a 
significant role and the rate is somewhat less sensitive to 
the details of the low-lying excitations. Therefore, even 
though the GT strength distribution calculated in QRPA 
is quite different from that observed in the experiment, 
the deduced EC rates coincidcntally match quite well. 



D. 



Ti 



Experimental information about GT transitions from 
the 51 V (| ) ground state to (§,|,§)~ final states in 
51 Ti is available from a 51 V(n,p) experiment [13] and a 
51 V(d, 2 He) experiment (48|. The excitation-energy reso- 
lution achieved in the (n,p) experiment was 0.9 MeV and 
the results are presented in Fig. E^a). Bins of 1-MeV wide 
were used, as in the original work. The results from the 
51 V(<i, 2 He) experiment are shown in Fig. [7jb). A reso- 
lution of 110 keV was achieved, allowing for a detailed 
extraction of B(GT) values up to E x rs 6.5 MeV. Here, 
we used Fig. 3 (top) in Ref. [48| to obtain B(GT) values 
on a state-by-state basis. 

Figs. EJc,d,e) show the theoretical GT strength dis- 
tribution from the QRPA (divided by a factor of 4) and 
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FIG. 5. GT transitions from 50 V to 50 Ti: (a) from (d, 2 He) 
data (a known but unobserved 6 + state at 3.198 MeV is indi- 
cated), (b) from calculations in QRPA (divided by a factor of 
3), (c) from shell- model calculations using the KB3G interac- 
tion and (d) from shell-model calculations using the GXPFla 
interaction. In (e), running sums of B(GT) as a function of 
excitation energy are plotted. 



shell-model calculations in full pf model space using the 
KB3G and GXPFla interactions, respectively. Fig. [TJf) 
shows the running sum of B(GT) as a function of exci- 
tation energy. 

The GT strength distributions calculated in the shell- 
model using the GXPFla and KB3G interactions agree 
well with those extracted from the data. Above 6.5 MeV, 
where (d, 2 He) strengths are not available, a considerable 
amount of strength is seen in the (n,p) data. The shell- 
model calculations predict some strength to be located at 
higher excitation energies, but slightly less than observed 
in the (n,p) data. The QRPA strength distribution ex- 
hibits little fragmentation and produces two strong GT 
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FIG. 6. Electron-capture rates for 50 V(e", ^ e ) 50 Ti based 
on experimental and theoretical GT strength shown in Fig. 
[5]G!T, plotted as a function of temperature at two densities, 
pY e =10 7 g/cm 3 (a) and pY e =10 9 g/cm 3 (b). 



transitions at ~ 2 and ~ 7 MeV. A weak transition is 
predicted at very low excitation energies, but not seen in 
the high-resolution 51 V(d, 2 He) data. The summed GT 
strength up to 11 MeV in the QRPA calculation is about 
40% higher than seen in the 51 Y(n,p) data. 

The close correspondence between the GT strengths 
calculated in the shell-model and those extracted from 
the 51 V(<i, 2 He) data results in a close match of the de- 
duced EC rates, as shown in Fig. [5] Since Qec = —2.982 
MeV and the bulk of the GT strength is located below 6 
MeV, the strength not extracted at high excitation ener- 
gies in the 51 V(d, 2 He) data is of little consequence, even 
at the higher density. Because of the limited resolution 
achieved in the 51 V(n,p) experiment, some strength is 
found at relatively low excitation energies, resulting in an 
increase of the EC rate, in particular at low temperatures. 
In the QRPA calculation, excess strength is predicted at 
low excitation energies. Combined with the prediction 
of a low-lying state not observed in the high-resolution 
51 V(<i, 2 He) data, the EC rates calculated based on the 
QRPA calculations are significantly higher than those 
based on that data, at both densities considered. 



E. 



1 Pe- 



4 Mn 



Information about GT transitions from the 54 Fe(0 + ) 
ground state to 1 + states in 54 Mn is available from a 
54 Fe(n,p) experiment, performed at 300 MeV [4!| [50| . 
The excitation-energy resolution was 1.2 MeV. A second 
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FIG. 7. GT transitions from 51 V to 51 Ti: (a) from (n,p) data, 
(b) from (d, 2 He) data, (c) from calculations in QRPA (divided 
by a factor of 4), (d) from shell- model calculations using the 
KB3G interaction and (e) from shell-model calculations using 
the GXPFIa interaction. In (f), running sums of B(GT) as a 
function of excitation energy are plotted. 



54 Fe(n,p) experiment was performed at 100 MeV [5lj . 
but since the resolution was relatively poor (2.8 MeV), 
we used the earlier data for the comparison with the the- 
oretical calculation. Those results are shown in Fig. ■ 
Since Table I in Ref. [Icj gives the B(GT) values in 2- 
MeV wide bins, we have instead used Fig. 10(a) and the 
quoted GT unit cross section of 5.1±0.8 mb/sr in that ref- 
erence to produce a GT strength distribution in 1-MeV 
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FIG. 8. Electron-capture rates for 51 V(e~, i/ e ) 51 Ti, based on 
experimental and theoretical GT strength distributions shown 
in Fig. [Jj plotted as a function of temperature at two den- 
sities, pY e =10 7 g/cm 3 (a) and pY e =10 9 g/cm 3 (b). Note 
that the rates based on the shell-model calculations with the 
KB3G and GXPFIa interactions overlap. 



wide bins. We note that a 54 Fe(<i, 2 He) experiment has 
also been performed [52| with slightly better resolution 
(~ 700 keV) than in the above-mentioned 54 Fe(n,p) ex- 
periments. However, the (n,p) data was used to extract 
information about GT excitations from the 54 Fe(<i, 2 He) 
experiment, which made it impossible to treat the lat- 
ter as an independent measurement. It is, therefore, not 
included in the present analysis. 

Figs. [HIb,c,d) show the theoretical calculations for the 
GT strength distribution from the QRPA calculation (di- 
vided by a factor of 5) and shell-model calculations in full 
pf model space using the KB3G and GXPFIa interac- 
tions, respectively. Fig. [Hie) shows the running sums of 
B(GT) as a function of excitation energy. 

Both shell-model calculations do well in describing the 
total GT strength up to an excitation energy of 10 MeV 
extracted from the 5i Fe(n,p) data, but the strength dis- 
tribution is spread out more in the calculation with the 
GXPFIa interaction than in the calculation with the 
KB3G interaction. In the 54 Fe(n,p) data a significant 
amount of strength is found between 0.5 MeV and 1.5 
MeV. In contrast, both shell-model calculations predict 
the lowest transitions to occur to 1 + states in 54 Mn at 
about 1.5 MeV, which is consistent with the location of 
positively identified 1 + states at 1.391 MeV and 1.454 
MeV [53|]. The results from the 54 Fe(d, 2 He) experiment 
[52j also seem to indicate that a large fraction of the 
strength found below 1.5 MeV in the (n,p) experiment 
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is in fact due to the transitions to the states located at 
E x k, 1.4 MeV. We conclude that the correspondence be- 
tween the shell-model calculations and the data would 
probably improve if high-resolution spectra were avail- 
able. The QRPA calculations do not reproduce the ex- 
perimental spectrum well, concentrating nearly all GT 
strength in the population of two states at E x « 4 MeV. 
The summed strength up to 10 MeV is 60% higher than 
observed in the data. Although not distinguishable in 
Fig. Efb), a weak (B(GT) = 0.016) transition to a state 
at E x = 0.55 MeV is also predicted. 

Given the relatively large amount of GT strength 
placed at low-excitation energies, the deduced EC rate 
from the 54 Fe(n,p) experiment is significantly higher 
than those deduced from the theoretical strength distri- 
butions, as shown in Fig. [TO] Because of the effects re- 
lated to the poor resolution of the data discussed above, 
this discrepancy is at least partially artificial. Due to the 
presence of the weak transition to the 1 + state at 0.55 
MeV predicted in the QRPA calculations, the EC rates at 
low density (Fig. HUTa)) are quite similar to those based 
on the shell-model calculations. However, the lack of any 
significant GT strength up to 3.7 MeV results in a large 
underestimation of the EC rate at the higher densities 
(Fig. flDT b)). in particular at low temperatures. 



'Mn^ 5j Cr 



The GT strength distribution for transitions from 
the 55 Mn(§~) ground state to (§,§,f)~ states in 55 Cr 
has been measured via the 55 Mn(n,p) reaction at 198 
MeV H|. Additionally, the B{GT) for the 55 Mn(f~, 

g.s.)— s- 55 Cr(| , g.s.) transition is known from 55 Cr(/3~) 
data [Hil [Iff. Fig. [TlTa) shows the experimental GT 
distribution, based on a combination of results from the 
55 Mn(n,p) and 55 Cr(/3") experiments. The B{GT) for 
the transition to the ground state was fixed to the mea- 
sured value (0.0242) from 55 Cr(/3"). Additional strength 
measured below 0.5 MeV in the 55 Mn(n,p) experiment 
was placed at E^( 55 Cr)=0.25 MeV. For the remainder 
of the spectrum, values from Ref. [HiJ were adopted. 
Figs. [TTTb-d) show the theoretical calculations for the 
GT strength distribution. Fig. II If el shows the running 
sums of B(GT) as a function of excitation energy. 

The shell-model calculations using the KB3G and 
GXPFIa interactions reproduce the experimental spec- 
trum moderately well, including the transition to the 
ground state of 55 Cr. In the QRPA calculations, too 
much strength is predicted below 2 MeV and the total 
strength up to 10 MeV is about 30% higher than in the 
data. In addition, no strength is predicted for the tran- 
sition to the 55 Cr ground state and the first transition 
populates a state at 0.51 MeV. 

The EC rates calculated by using the experimental GT 
strength distribution is slightly higher than the EC rates 
based on the shell-model calculations (see Fig. IT2"1) . How- 
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FIG. 9. GT transitions from 54 Fe to 54 Mn: (a) from (n,p) 
data, (b) from calculations in QRPA (divided by a factor of 
5), (c) from shell- model calculations using the KB3G interac- 
tion and (d) from shell-model calculations using the GXPFIa 
interaction. In (e), running sums of B(GT) as a function of 
excitation energy are plotted. 



ever, it is difficult to draw very strong conclusions from 
these minor discrepancies because of the limited resolu- 
tion of the (n,p) experiment. Since Qec = —3.114 MeV, 
the calculated EC rates are very sensitive to the strength 
distributions at low excitation energies, in particular at 
the lower density (Fig. IT27 a)). The EC rates based on 
the QRPA calculations are slightly higher than the rates 
based on the data and the shell-model calculations (ex- 
cept at the lower density and the lowest temperatures), 
due to the presence of the two strong transitions below 2 
MeV, which are not observed experimentally and do not 
appear in the shell-model calculations. 
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FIG. 10. Electron-capture rates for 54 Fe(e~, v e ) 54 Mn based 
on experimental and theoretical GT strengths shown in Fig. 
[9] plotted as a function of temperature at two densities, 
pY e =10 7 g/cm 3 (a) and pY e =10 9 g/cm 3 (b). 



G. 



J Fe- 



3 Co 



Information about GT transitions from the ground 
state of 56 Fe(0 + ) to 56 Mn is available from a 56 Fe(n,p) 
experiment, performed at 198 MeV The excitation- 
energy resolution was 1.3 MeV. A second 56 Fe(n,p) ex- 
periment performed at 100 MeV [5l[ had a relatively poor 
excitation-energy resolution of 2.8 MeV. Hence, we used 
the earlier data which is shown in Fig. fTBTa). In the cur- 
rent analysis, strength reported below E x = MeV in 
Ref. [54[ was added to the first energy bin. It should be 
noted that high-resolution (110 keV FWHM) 56 Fe(d, 2 He) 
data are available [57j . but GT strengths were not pre- 
sented in that reference. Figs. ITBT b-d') show the theoret- 
ical calculations for the GT strength distribution. Fig. 
ITBT e) shows the running sums of B{GT) as a function of 
excitation energy. 

Taking into acount the limited experimental resolution 
of the 56 Fe(n,p) data, the shell-model calculations match 
well with the data up to an excitation energy of about 
7 MeV. Above that energy, significant amounts of GT 
strength were extracted from the (n,p) data, whereas no 
strength is present in the calculated distributions. The 
QRPA calculations do not reproduce the experimental 
strength distribution well, missing the strength below 
2 MeV and overestimating the total strength by about 
50%. Although GT strengths are not available from the 
high resolution 56 Fe(d, 2 He) experiment (57|, that data set 
clearly shows two strong transitions to 1 + states at 0.11 
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FIG. 11. GT transitions from 55 Mn to 55 Cr: (a) from (n,p) 
and and /3— decay data, (b) from calculations in QRPA, (c) 
from shell-model calculations using the KB3G interaction and 
(d) from shell-model calculations using the GXPFIa interac- 
tion. In (e), running sums of B(GT) as a function of excita- 
tion energy are plotted. 



MeV and 1.2 MeV. The shell- model calculations with 
both interactions exhibit a similar feature, although the 
calculation with the KB3G interaction locates the first 
1 + state slightly too high in excitation energy. 

The EC rates based on the GT strengths calculated 
in the shell-models fall slightly below those based on the 
strengths extracted from the 56 Fe(n,p) as shown in Fig. 
1141 The minor discrepancies are likely enhanced by the 
poor resolution of the {n,p) data, which results in the 
placement of some GT strength at artificially low exci- 
tation energies. Since the excitation energy of the first 
state is slightly higher in the calculation with the KB3G 
interaction than with the GXPFIa interaction, the EC 
rates are slightly lower. Since the excitation energy of the 
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FIG. 12. Electron-capture rates for 55 Mn(e", ^ e ) 55 Cr based 
on experimental and theoretical GT strengths shown in Fig. 
1111 plotted as a function of temperature at two densities, 
pY e = W 7 g/cm 3 (a) and pY e =10 9 g/cm 3 (b). 

first state that can be captured into is much higher in the 
QRPA calculations than in the data or the shell-model 
calculations, the associated EC rates are much lower. At 
increasing temperatures, when the strengths at higher 
excitation energies play a more prominent role, the dis- 
crepancies reduce, assisted by the fact that the total GT 
strength predicted in the QRPA calculation is about 50% 
higher than in the data or the shell-model calculations. 

H. 58 Ni^ 58 Co 

Experimental information on GT transitions from the 
58 Ni(0 + ) ground state to 1 + states in 58 Co is available 
from 58 Ni(n,p) Q, 58 Ni(d, 2 He) \&&\$ and 58 Ni(£, 3 He) 
[60] | experiments. The results are displayed in Figs. 
[PoTa). (b) and (c), respectively. The energy resolu- 
tions in the 58 Ni(n,p), 58 Ni(d, 2 He) and 58 Ni(i, 3 He) data 
sets were 1.2 MeV, 130 keV and 250 keV, respectively. 
The 58 Ni(n,j?) data was available in energy bins of 1 
MeV. In the current analysis, strength recorded below 
E x — MeV was added to the first energy bin. For 
excitation energies below 4 MeV, GT strengths from the 
58 Ni(d, 2 He) experiment were available on a state-by-state 
basis. Above that excitation energy, an energy binning of 
1 MeV was used. The GT strengths from the 58 Ni(i, 3 He) 
experiment were available in 250-keV wide bins up to an 
excitation energy of 10 MeV. Figs. rToT d-f) show the theo- 
retical calculations for the GT strength distribution from 
the QRPA (divided by a factor of 6) and shell-model cal- 



data Fe(n,p) 




E x ( Mn) (MeV) 

FIG. 13. GT transitions from 56 Fe to 56 Mn: (a) from (n,p) 
data, (b) from calculations in QRPA (divided by a factor of 
3), (c) from shell- model calculations using the KB3G interac- 
tion and (d) from shell-model calculations using the GXPFIa 
interaction. In (e), running sums of B(GT) as a function of 
excitation energy are plotted. 



culations in truncated pf model space (5 holes in the 
shell) using the KB3G and GXPFIa interactions, respec- 
tively. Fig. 115( g) shows the running sums of B(GT) as 
a function of excitation energy. 

The main difference between the shell-model calcula- 
tions and the experimental data occurs at an excitation 
energy of about 2 MeV. In the calculation with the KB3G 
interaction, too much strength is assigned to a single 
transition, whereas in the calculation with the GXPFIa 
interaction, too little strength is present near that ex- 
citation energy. Since slightly more strength is present 
at higher excitation energies in the calculation with the 
GXPFIa interaction, the summed strengths up to 10 
MeV for the two shell-model calculations are very close 
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FIG. 14. Electron-capture rates for 56 Fe(e~, v e ) 56 Mn based 
on experimental and theoretical GT strengths plotted in Fig. 
1131 as a function of temperature at two densities, pY e =10 7 
g/cm 3 (a) and pY e = 10 9 g/cm 3 (b). 



and also consistent with the experimental results. In the 
QRPA calculations, most of the strength is concentrated 
around an excitation energy of 4 MeV and the summed 
strength is a factor of two higher than observed in the 
data. A very weak transition (B(GT) = 0.00321) pop- 
ulates a 1 + state at 0.21 MeV; it is not visible in Fig. 
HSId). 

The EC rates calculated from the three data sets match 
each other quite well, as shown in Fig. Q1)J Only at the 
lower density and at low temperatures (Fig. ITST a)) is the 
rate based on the 58 Ni(n,p) data significantly higher than 
those based on the 58 Ni(d, 2 He) and 58 Ni(£, 3 He) data sets. 
This is likely due to the relatively poor energy resolution 
achieved in the 58 Ni(n,p) experiment, which results in 
the placement of a small amount of strength at excitation 
energies below 1 MeV. 

At ( oY e =10 7 g/cm 3 , the EC rates based on GT 
strengths predicted in the shell-model calculations do 
well in reproducing those based on the GT strengths ex- 
tracted from the three data sets. At the higher den- 
sity, the EC rate deduced from the calculation with the 
KB3G interaction is too high, due to the fact that too 
much strength is placed at relatively low excitation ener- 
gies. For the EC rate deduced from the calculation with 
the GXPFIa interaction the reverse is the case, although 
the discrepancy with the rates based on the experimental 
data is smaller and vanishes at the highest temperatures. 
EC rates based on the QRPA calculations are too low, 
except for the low-temperature region at the lower den- 
sity (due to the predicted weak transition to a state at 



0.21 MeV) and for the high-temperature region at the 
higher density (where the strong states near 4 MeV start 
playing a large role). 



I. 59 Ctw 59 Fe 

The GT strength distribution for transitions from the 
59 Co(| ) ground state to (§,|,§)~ states in 59 Fe has 
been measured via the 59 Co(n,p) reaction at 198 MeV 
(47| . The excitation-energy resolution was 0.9 MeV 
and the data were presented in 1-MeV wide bins. GT 
strength recorded below E x = MeV were added to the 
0-1 MeV bin for the present study. The experimental 
strength distribution is shown in Fig. UTT a). Figs. [TWb- 
d) show the theoretical calculations for the GT strength 
distribution. Fig. I17f e) shows the running sums of 
B(GT) as a function of excitation energy. 

A reasonable correspondence between the data and the 
theoretical calculations using the shell models is found. 
The calculation with the KB3G interaction has a rela- 
tively strong transition to a state at 2.4 MeV, which is 
weaker in the calculations with the GXPFIa interaction. 
In both shell-model calculations the summed strength up 
to an excitation energy of 6 MeV is about 50% higher 
than observed in the data. The discrepancy in the 
summed strength becomes smaller at higher excitation 
energies, because of the relatively large amount of GT 
strength found above 6 MeV in the data. In the QRPA 
calculation a strong transition to a state at 1.5 MeV is 
predicted, but it is not present in the data. The remain- 
der of the strength appears at around 6 MeV, about 2 
MeV higher than where the bulk of the GT strength was 
found in the experiment. The summed strength up to 
10 MeV is overestimated by about 60% in the QRPA 
calculation. 

At pY e =10 7 g/cm 3 , and at temperatures below 5 x 10 9 
K, the EC rates calculated based on the experimentally 
extracted GT strength distribution are higher than those 
deduced from the theoretical strength distributions (see 
Fig. [T8Fa)). This could, in part, be due to the placement 
of some strength at artificially low excitation energies in 
the experimental distribution as a consequence of the lim- 
ited resolution. Even the strong transition to the state at 
1.5 MeV in the QRPA calculations does not play a large 
role at these low densities and temperatures. However, 
at /3Y e =10 9 g/cm 3 (see Fig. [T^b)), this transition dom- 
inates the EC rate, which then exceeds those based on 
the data and the shell-model calculations. With slightly 
more strength placed at lower excitation energies, the EC 
rates based on the shell-model calculation with the KB3G 
interaction provides a better match to the EC rates de- 
duced from the experimental strength distribution than 
the EC rates based on the shell-model calculation with 
the GXPFIa interaction, except at the highest stellar 
temperatures. 
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FIG. 15. GT transitions from 58 Ni to 58 Co: (a) from (n,p) 
data, (b) from (d, 2 He) data, (c) from (t, 3 He) data (d) from 
calculations in QRPA (divided by a factor of 6), (e) from 
shell-model calculations using the KB3G interaction and (f) 
from shell-model calculations using the GXPFIa interaction. 
In (g), running sums of B(GT) as a function of excitation 
energy are plotted. 
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FIG. 16. Electron-capture rates for 58 Ni(e~ , i/ e ) 58 Co based on 
experimental and theoretical GT strengths from Fig. 1151 plot- 
ted as a function of temperature at two densities, pY e — 10 7 
g/cm 3 (a) and pY e =10 9 g/cm 3 (b). 



The GT strength distribution for transitions from the 
60 Ni(0 + ) ground state to 1 + states in 60 Co has been mea- 
sured via the 60 Ni(n,p) reaction at 198 MeV [gl|. The 
energy resolution was 860 keV and the strength distri- 
bution was presented in 1-MeV wide excitation-energy 
bins. The spectrum below E x ( 60 Co) = 3.5 MeV was 
reanalyzed in Ref. (62j and in the present work the 3- 
Gaussian fit presented in Fig. 6(b) of Ref. [62[ was 
used to represent the data. For E x ( &0 Co) > 3.5 MeV, 
the extracted strengths from Ref. [6l| were used. The 
strength distribution is shown in Fig. HUT a). Further 
experimental information is available from a 60 Ni(p, n) 
experiment at 134 MeV [62], in which GT transitions to 
states in 60 Cu at high excitation energies were studied. 
These states are the T = 3 analogs of the low- lying 1 + 
states in 60 Co. Here, we used the extracted B(GT) val- 
ues from a 3-Gaussian fit (see Fig. 6(a) of Ref. (62[) 
as shown in Fig. rT97 b). Since the T — 3, 1 + states in 
60 Cu are situated in a continuum of T = 1 and T = 2 
states, extraction of higher-lying GT strengths from the 
60 Ni(p, n) experiment was not possible. 

Figs. [TWc.d.e) show the theoretical calculations for the 
GT strength distribution from the QRPA (divided by a 
factor of 2) and shell-model calculations in truncated pf 
model space (5 holes in the ff/ 2 shell) using the KB3G 
and GXPFIa interactions, respectively. Fig. IT97 shows 
the running sums of B(GT) as a function of excitation 
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FIG. 17. GT transitions from 59 Co to 59 Fe: (a) from (n,p) 
data, (b) from calculations in QRPA (without additional 
scaling), (c) from shell-model calculations using the KB3G 
interaction and (d) from shell-model calculations using the 
GXPFIa interaction. In (e), running sums of B(GT) as a 
function of excitation energy are plotted. 



energy. 

Discrepancies observed between the two available data 
sets for the transitions to states below 3 MeV make it 
difficult to draw strong conclusions on the quality of the 
theoretical calculations. Nevertheless, several observa- 
tion can be made. In the shell-model calculations with 
the KB3G interaction, a strong transition to a state at an 
excitation energy of 1.3 MeV is predicted, with a strength 
comparable to that of the transition to the state observed 
at 0.7 MeV in both data sets. In the calculation with 
the GXPFIa interaction, the excitation energy of this 
state is very close to the experimental value, but the 
strength lower by about a factor of 2. Several weaker 
transitions appear to states at excitation energies just 
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FIG. 18. Electron-capture rates for 59 Co(e~ , i/ e ) 59 Fe based 
on experimental and theoretical GT strengths plotted in Fig. 
1171 as a function of temperature at two densities, pY e — 10 7 
g/cm 3 (a) and pY e = 10 9 g/cm 3 (b). 



above the strong transition in the calculation with the 
GXPFIa interaction. The summed GT strengths up to 
an excitation energy of 8 MeV in the two shell-model cal- 
culations are consistent with that found in the 60 Ni(n,p) 
data. In the QRPA calculations almost all strength is 
associated with transitions to states with excitation en- 
ergies between 2 and 3 MeV and the summed strength 
exceeds that extracted from the data by about a factor 
of two. 

Because Q E c = -3.334 MeV. the EC rates on 60 Ni arc 
dominated by the GT transitions to the final states at low 
excitation energies. Therefore, in spite of the fact that no 
strength could be extracted from the 60 Ni(p, n) data at 
high excitation energies, the difference between the EC 
rates calculated based on the 60 Ni(n,p) and 60 Ni(j>, n) 
data sets are small, as shown in Fig. [501 Since the loca- 
tion of the final state at the lowest excitation energy is 
best reproduced in the shell-model calculation with the 
GXPFIa interaction, the deduced EC rates are also close 
to those deduced from the data. Although partially mit- 
igated by the overestimated strength of the transition to 
that state, the EC rates based on the GT strength dis- 
tribution calculated in the shell-model with the KB3G 
interaction is lower because the state is placed at higher 
excitation energies. For the rates based on the QRPA 
calculations, that effect is enhanced since most of the 
strength is located at even higher excitation energies. At 
higher densities and temperatures, where the rates be- 
come more sensitive to the total GT strength present in 
the spectra, the experimental and theoretical rate esti- 
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FIG. 19. GT transitions from 60 Ni to 60 Co: (a) from (n,p) 
data, (b) from T = 3 states observed in (p,n) data, (c) from 
calculations in QRPA (divided by a factor of 2), (d) from 
shell-model calculations using the KB3G interaction and (e) 
from shell-model calculations using the GXPFIa interaction. 
In (f), running sums of B(GT) as a function of excitation 
energy are plotted. 
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FIG. 20. Electron-capture rates for 60 Ni(e" , v e f°Co for ex- 
perimental and theoretical GT strengths plotted in Fig. I19I 
as a function of temperature at two densities, pY e =10 7 g/cm 3 
(a) and pY e =10 9 g/cm 3 (b). 



Experimental data available for GT transitions from 
the 62 Ni(0 + ) ground state to 1 + states in 62 Co are simi- 
lar in nature to those available for 60 Ni. 62 Ni(n,p) data 
at 198 MeV [gU provides the GT strength distribution 
extracted via a multipole decomposition analysis in bins 
of 1 MeV. Ref. [62] gives a reanalysis of the low-lying part 
(E x ( 62 Co) < 2.5 MeV) of the 62 Ni(n, p) excitation-energy 
spectrum (see Fig. 5(d) of Ref. (62[) and the deduced 
B(GT) values were used to generate the GT strength dis- 
tribution shown in Fig. HUa). For £ x ( 62 Co) > 2.5, the 
GT strengths extracted from the multipole decomposi- 
tion analysis in Ref. [6l| were used. Further experimental 
information is available from a 62 Ni(p, n) experiment at 



mates converge. 



134 MeV [62j. GT transitions to states in 62 Cu at high 
excitation energies were studied. These states are the 
T = 4 analogs of the low-lying 1 + states in 62 Co. Here, 
we used the extracted B(GT) values from a 2-Gaussian 
fit (see Table II of Ref. [Hj]) as shown in Fig. (Hfb). 
Since the T = 4, 1 + states in 62 Cu are situated in a con- 
tinuum of T — 2 and T = 3 states, extraction of higher- 
lying GT strengths from the 62 Ni(p, n) experiment was 
not possible. Figs. [2H c-e) show the theoretical calcula- 
tions for the GT strength distribution. Fig. l2"lT f) shows 
the running sums of B(GT) as a function of excitation 
energy. 

The experimental results from the two sources of data 
are consistent in the excitation-energy range where both 
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are available. Relatively little GT strength was reported 
at excitation energies above 3 MeV in the analysis of the 
62 Ni(n,p) data. Both shell- model calculations correctly 
predict a strong transition to a final state at low excita- 
tion energies. The predicted excitation energy in the cal- 
culation with the KB3G interaction is too high by about 
0.4 MeV and 50% too strong compared to the experi- 
mental strength. Although a little too weak (by about 
20%), the calculation with the GXPFla interaction pre- 
dicts the excitation energy correctly. Both shell-model 
calculations predict little strength at higher excitation 
energies and in the case of the KB3G calculation some 
of the excess strengths predicted for the transition to 
the first 1 + state is balanced by the reduced amount of 
strength predicted at higher excitation energies. Conse- 
quently, in both calculations the total GT strength is in 
close correspondence with the data. In the QRPA calcu- 
lation, a strong state at low excitation energies is absent 
and the bulk of the strength, in excess of the data by 
about a factor of 2, is concentrated in a narrow excita- 
tion energy window centered around 2 MeV. There are, 
however, about 10 very weak transitions to states below 
1.0 MeV (with excitation energies of as low as 0.06 MeV), 
which are not all visible in Fig. [2TTc'l . 

Since Qec — —5.862 MeV, the EC rates are strongly 
dominated by transitions to final states at low excitation 
energies, even at the higher density of 10 9 g/cm 3 . For 
the strength distributions extracted from the data and 
the strength distributions calculated in the shell-models, 
the response is dominated by the strong transition to the 
first 1 + state. The slightly higher excitation energy of 
this state in the calculation with the KB3G interaction 
reduces the deduced EC rate compared to that deduced 
from the calculation with the GXPFla interaction. How- 
ever, this reduction is countered by the slightly higher 
strength for the transition to this state in the calculation 
with the KB3G interaction. 

The opposite effect is seen for the EC rate based on 
the QRPA calculation due to several very weak transi- 
tions to states below 1 MeV. The effects of the small 
transition strengths on the EC rates are countered by 
the low-excitation energies of the final states and the de- 
duced EC rates are relatively close to the other EC rates, 
as shown in Fig. 
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Data on GT transitions from 64 Ni are available from 
64 Ni(n,p) [61| and 64 Ni(d, 2 He) [63| experiments. In ad- 
dition, the B(GT) for the transition from the 64 Ni(0+) 
ground state to the 64 Co(l + ) ground state can be de- 
duced from a 64 Co(£~) experiment (B(GT) = 0.627) 
(64| . The resolution in the 64 Ni(n,p) experiment was 
800 keV and the GT strength distribution was extracted 
through a multipole decomposition analysis in 1-MeV 
wide excitation-energy bins. In the present study, these 
data were combined with the GT strength extracted for 
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FIG. 21. GT transitions from 62 Ni to 62 Co: (a) from (n,p) 
data, (b) from T = 4 states observed in (p,n) data, (c) from 
calculations in QRPA (divided by a factor of 2), (d) from 
shell-model calculations using the KB3G interaction and (e) 
from shell-model calculations using the GXPFla interaction. 
In (f), running sums of B(GT) as a function of excitation 
energy are plotted. 



the 64 Co(/3~) transition to the ground state, as shown 
in Fig. I2"37 a). The energy resolution achieved in the 
64 Ni(<i, 2 He) experiment was 110 keV, allowing for extrac- 
tion of the B(GT) values on a state-by-state basis up to 
an excitation energy in 64 Co of 1.5 MeV and in narrow 
energy-bins up to 3.7 MeV (Fig. [21b)). The 64 Co(/?-) 
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FIG. 22. Electron-capture rates for 62 Ni(e~, i/ e ) 62 Co for ex- 
perimental and theoretical GT strengths, plotted as a func- 
tion of temperature at two densities, pY e = 10 7 g/cm 3 (a) and 
pY e = 10 9 g/cm 3 (b). 



data was used in Ref. [63[ to calibrate the extraction of 
GT strengths. The B(GT) for the transition from the 
64 Ni(0 + ) ground state to the 64 Co(l + ) ground state is 
therefore by definition equal to that deduced from the 
/3~ experiment [64]]. 

Figs. l23T c-e) show the theoretical calculations for the 
GT strength distribution. Fig. H2Jf) shows the running 
sums of B(GT) as a function of excitation energy. 

The dominant GT transition to the ground state is well 
reproduced (difference of less than 10%) by the shell- 
model calculation using the GXPFla interaction and 
overestimated (by about 40%) in the calculation using 
the KB3G interaction. Both shell-model calculations pre- 
dict less strength at higher excitation energies than ob- 
served in the data. The results from the high-resolution 
(d, 2 He) experiment suggest that the extraction of GT 
strengths from the (n,p) data using a MDA led to the 
loss of a small fraction of GT strength, at least up to 3.7 
MeV. Since GT strengths were not extracted from the 
(<i, 2 He) data above that excitation energy, is not clear 
whether that trend would continue up to the higher ex- 
citation energies. Even under the assumption that no 
strength was lost in the analysis of the (n,p) data, the 
total GT strength predicted in the shell-model calcula- 
tions is about 30% lower than observed in the data. The 
QRPA calculations fail to reproduce the strong transition 
to the ground state of 64 Co; most of the GT strength 
is located at excitation energies between 1 and 2 MeV 
and the total strength is in excess of the experimental 
strength by about a factor of two. 



Since Qec = -7.818 MeV, the EC rates on 64 Ni 
are very small and dominated by the transition to the 
ground state. Since the shell-model calculation with the 
GXPFla interaction provides a good match to the data, 
the EC rates deduced from this model are very close to 
those based on the experimental GT strength distribu- 
tion (see Fig. |2"I| . The EC rates based on the shell- model 
calculation with the KB3G interaction is too high by 
about 40% due to the overestimation of the GT strength 
of the transition to the ground state of 64 Co; that dis- 
crepancy is not visible due to the logarithmic scales of 
Fig. [M] EC rates based on the GT strength distribu- 
tion calculated in QRPA are too low, in particular at low 
temperatures, due to the absence of the strong transition 
to the ground state. 
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Data on GT transitions from the 64 Zn(0 + ) ground 
state to 64 Cu(l + ) states are available from 64 Zn(<i, 2 He) 
[65| and 64 Zn(i, 3 He) [HJ experiments. Additionally, 
the B(GT) for the transition from the 64 Zn(0 + ) ground 
state to the 64 Cu(l + ) ground state can be deduced from 
64 Cu(/3") experiments (B(GT) = 0.059) 66]. The en- 
ergy resolution achieved in the 64 Zn(d, 2 He) experiment 
was 115 keV, allowing for extraction of the B(GT)s on 
a state-by-state basis up to an excitation energy in 64 Cu 
of 5 MeV, as shown in Fig. l25Ta). The 64 Cu(/3") data 
was used in Ref. (6f| to calibrate the extraction of GT 
strengths and the B(GT) for the transition from the 
64 Zn(0 + ) ground state to the 64 Cu(l + ) ground state was, 
therefore, by definition equal to that deduced from the 
/3~ experiment. The energy resolution achieved in the 
64 Zn(i, 3 He) experiment was 280 keV and a multipolc 
decomposition analysis was performed for 250-keV wide 
excitation-energy bins up to 9 MeV. The extraction of 
the GT strengths in the analysis did not rely on a cali- 
bration of the unit cross section using the 64 Cu(/3~) data. 
However, the B(GT) for the ground-state to ground- 
state transition was set equal to that obtained from the 
64 Cu(^-) data, as shown in Fig. 125Tb). Figs. H3Jc-e) 
show the theoretical calculations for the GT strength dis- 
tribution. Fig. [25Tf) shows the running sums of B(GT) 
as a function of excitation energy. 

Compared to the other nuclei studied in this work, the 
GT strength distributions calculated in the shell-models 
provide a poor description of the data. The calcula- 
tion with the KB3G interaction places the bulk of the 
strength at much too high excitation energies and the 
calculation with the GXPFla interaction also places too 
little strength at low excitation energies. The strength 
calculation based on the QRPA calculations suffers from 
a similar problem. In addition, it overestimates the total 
strength by nearly a factor two. Both shell-model calcu- 
lations underestimate the GT strength for the transition 
to the ground state; by 50% for the GXPFla and and 
by 80% for the KB3G interaction. The GT strength for 
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FIG. 23. GT transitions from 64 Ni to 64 Co: (a) from (n,p) 
and /9-decay data, (b) from (d, 2 He) and /3-decay data, (c) from 
calculations in QRPA, (d) from shell-model calculations using 
the KB3G interaction and (e) from shell-model calculations 
using the GXPFIa interaction. In (f), running sums of B(GT) 
as a function of excitation energy are plotted. 



the transition to the ground state is underestimated by 
a factor of 16 in the QRPA calculations and not visible 
in the plotted spectrum. 

Not surprisingly, the EC rates based on the theoreti- 
cal calculations underestimate those based on the experi- 
mental strength distributions. Since Qec = —1-09 MeV, 
GT strengths at relatively high excitations energies can 
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FIG. 24. Electron-capture rates for 64 Ni(e~ , i/ e ) 64 Co for ex- 
perimental and theoretical GT strengths plotted in Fig. I23I 
as a function of temperature at two densities, pY e =10 7 g/cm 3 
(a) and pY e =10 9 g/cm 3 (b). 



contribute significantly, especially at the higher density. 
This somewhat reduces the discrepancies between the EC 
rates based on the data and the theory. 



IV. DISCUSSION 

Based on the discussion of the GT strength distribu- 
tion for individual nuclei in the previous section, one 
can conclude that the shell-model calculations with the 
KB3G and GXPFIa interaction describe the main fea- 
tures of experimental GT strength distributions. The 
exceptions are the cases of 64 Zn and 45 Sc. For the latter 
case, high-resolution, background-free data are needed 
to better assess the discrepancies between experiment 
and shell-model theory. For nuclei with mass number 
A > 56, the GT strength distributions calculated with 
the GXPFIa interaction exhibit more fragmentation and 
are somewhat broader than those calculated with the 
KB3G interaction. With the latter interaction, more 
strength is placed at low-excitation energies, and less at 
higher excitation energies. Similar differences between 
these two shell-model calculations also became apparent 
in the recent measurement of the GT strength distribu- 
tion from unstable 56 Ni [67|. In general, the calculations 
performed in the QRPA formalism do not describe the 
details of the strength distributions well, and the total 
strength is overestimated by up to a factor of two. 

To gain more quantitative insight in the impact of the 
differences between the experimental and theoretical GT 



21 



0.5 



0.5 



0.5 



0.5 



0.5 



data 64 Zn(d, 2 He)+ 64 Cu(P") 



(a) 



II I 

data 64 Zn(t, 3 He)+ 64 Cu(p") 



(b) 



il I . I I i 

QRPA/1 



I h .I 



Hill I.I 



KB3G 



(d) 



ill i . I i II 



GXPFIa 



data^ 4 Zn(d 2 He) + a4 Cu(p-) 
data 64 Zn(t, He)+ 64 Cu(P") 



(f) 



KB3G 
GXPFIa 
QRPA/1 



•Jri. 



E x ( 64 Cu) (MeV) 

FIG. 25. GT transitions from 64 Zn to 64 Cu: (a) from (d, 2 He) 
and /3-decay data, (b) from (i, 3 He) and /3-decay data, (c) from 
calculations in QRPA, (d) from shell-model calculations using 
the KB3G interaction and (e) from shell-model calculations 
using the GXPFIa interaction. In (f), running sums of B(GT) 
as a function of excitation energy are plotted. 



strength distributions on the derived EC rates, the EC 
rates on the ground state for each of the 13 p/-shell nuclei 
investigated in this work are compared in Fig. [57] for two 
specific stellar density-temperature combinations (cases 
I and II: see section|TT|. In Fig. l2"Tf a) rates at a density of 
pY e —l{) 7 g/cm 3 and a temperature of T = 3xl0 9 K are 
compared (case I), and in Fig. l2"Wb). rates at a density 
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FIG. 26. Electron-capture rates for 64 Zn(e", !/ e ) 64 Cu for 
experimental and theoretical GT strengths shown in Fig. 
I25I plotted as a function of temperature at two densities, 
pY e =10 7 g/cm 3 (a) and pY e =10 9 g/cm 3 (b). 



of pY e =10 9 g/cm 3 and a temperature of T = 10xl0 9 
K are compared (case II). Depending on the value of 
Qr = Ei — Qec, the magnitudes of the absolute EC rates 
varies strongly between the cases studied here. There- 
fore, we chose to plot relative EC rates, and used the 
rates derived from the theoretical strength distributions 
using the shell-model with the KB3G interaction as a 
reference. 

For case I (Fig. [2"7T a)). the EC rates are very sensi- 
tive to the details of the GT strength distribution at low 
excitation energies. Of the nuclei studied, only in the 
cases of 45 Sc and 64 Zn is e F (T = 0) = 1.2 MeV> Q R 
for the transition to the states at the lowest excitation 
energies. In all other cases EC can only proceed at finite 
temperature. The placement of the states thus becomes 
critically important because the phase space for EC de- 
creases rapidly with increasing excitation energy. 

Since the resolution of the GT strength distributions 
extracted from the (n,p) experiments is limited, strength 
might be placed at artificially low excitation energies, re- 
sulting in an overestimate of the EC rates. This is clearly 
visible in Fig. |27Ja) for the cases of 48 Ti, 51 V and 58 Ni, 
for which (n,p) and high-resolution data are available. 
For the EC rates derived from GT strengths extracted 
from the 60 ' 62 ' 64 Ni(n,p) data, such artificial inflation does 
not occur as the GT strength at low excitation energies 
is completely dominated by a transition to a single state. 
Even with poor experimental resolution, the excitation 
energy of this state could be established with little un- 
certainty. 
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Ignoring the cases for which the excitation energies of 
the states at low excitation energies could not be well 
determined, it is clear that the deviations between the 
EC rates based on the shell-model calculations and the 
EC rates based on experimental GT strength distribu- 
tions are relatively small compared to the deviations ob- 
served for the rates based on the QRPA calculations. 
To quantify the differences, we defined an average devi- 
ation (Aec) between the EC rates based on the charge- 
exchange data (X exp ) and the theory (X th ), as a fraction 
of the experiment values as: 

^ N \th _ y^xp 

Aec = nY1 1 x expl ' (7) 

i— 1 1 

where the sum refers to all nuclei for which high resolu- 
tion data are available, or the GT strength distribution 
at low excitation energy is otherwise well-defined: 48 Ti, 
51 V, 58 Ni, 60 Ni, 62 Ni, 64 Ni and 64 Zn were used. Although 
high-resolution data are available for 50 V, we excluded 
them from this analysis because of the ambiguity related 
to the strength of the GT transition to the first 6 + state 
in the experimental analysis. 

Since Aec can be small when averaged over several 
nuclei if pos itive a nd negative deviations cancel out, we 
also defined |Asc|: 

\A E c \ = -y 1 ^exy > (8) 

which represents the average value of the absolute devia- 
tion between the EC rates based on the charge-exchange 
data (using the same 7 nuclei as for the calculation of 
Aec) and the theory, as a fraction of the EC rate based 
on the data. The values of Aec and |A.ec| are presented 
in Table HTl for each of the three theoretical models tested 
in this work, and for both temperature-density combina- 
tions defined above. 

For case I, the rates based on the shell-model calcula- 
tions with the GXPFla and KB3G interactions under- 
estimate the rates based on the data on average by ap- 
proximately 30% (Aec ~ 0.3). The average absolute 
deviation (|A_bc|) is slightly lower for the calculation us- 
ing the GXPFla interaction compared to the calculation 
using the KB3G interaction. However, given the uncer- 
tainties in the experimental strengths and rates, one can- 
not draw very strong conclusions based on this difference. 
The rates based on GT strength distributions produced 
in QRPA deviate strongly from those based on the data 
and it is clear that at low stellar densities and tempera- 
tures, the QRPA rates should be used with caution when 
modeling late stellar evolution. 

For case II, derived EC rates are less sensitive to the 
details of the strength distribution and the differences 
between the rates calculated based on theoretical and ex- 
perimental strength distributions become smaller, as can 
be seen in Fig. l2"7T b) (notice the difference in the scale 
of the axis depicting the rates from Fig. l2"7T a)). The EC 



rates derived from (n, p) data, for which the GT strength 
at low-excitation energies could not be accurately placed, 
are still too high compared to those derived from the 
high resolution data. However, the differences have re- 
duced to less than a factor of two. The EC rates de- 
rived from the shell-model calculations both match well 
with those derived from the data. The average deviation 
(Aec) between the EC rates derived from these calcula- 
tions and the data for the 7 nuclei, for which the strength 
at low excitation energies could be well-determined from 
the experiments, is less than 10% (see Table ITT|) . For the 
calculations with th e GXP Fla interaction, the average 
absolute deviation (lA^cl) is also below 10%. For the 
strength extracted from the calculations with the KB3G 
interaction the average absolute deviation is 30%. 

Since the EC rates at higher densities and tempera- 
tures are more sensitive to the overall GT strength, rather 
than the exact location of the daughter states at low exci- 
tation energies, the EC rates calculated based on the GT 
strength distributions from the QRPA formalism deviate 
much less from the rates based on the charge- exchange 
data than at the lower density and temperature. On av- 
erage, the rates are about 70% too high, which roughly 
corresponds with the average overestimation of the total 
GT strength in this theory, in comparison to the data. 



V. SUMMARY AND CONCLUSION 

We have evaluated 25 sets of data available from (n,p), 
(d, 2 He), (<, 3 He) and (p,n) experiments, providing infor- 
mation about GT strength distributions in the AT Z = +1 
direction for 13 nuclei in the pf shell. Three sets of the- 
oretical strength distributions (two of which were cal- 
culated in shell-models using the KB3G and GXPFla 
interactions and one was based on QRPA calculations 
using ground-state deformation parameters and masses 
from the finite-range droplet model) were tested against 
the experimental results. In addition, EC rates of rel- 
evance to late stellar evolution were derived from the 
experimental and theoretical GT strength distributions 
and compared as well. 

The GT strength distributions calculated in the shell- 
models do about equally well in reproducing the data, 
whereas large differences were observed for the QRPA 
calculations. Consequently, derived EC rates from the 
shell-model calculations are also much closer to the EC 
rates derived from the experimental GT strength dis- 
tributions than those calculated on basis of the QRPA 
framework. In order to perform quantitative comparisons 
between experimental and theoretical EC rates, we only 
used nuclei for which high-resolution data exist or for 
which the location of daughter states at low excitation en- 
ergies could be established with small uncertainty. At rel- 
atively low stellar densities and temperatures (pY e =10 
g/cm 3 and T = 3xl0 9 K), the EC rates derived from GT 
strengths calculated in the shell models were on average 
about 30% higher than the EC rates derived from the 
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FIG. 27. (color online) Comparison of EC rates for 13 p/-shell nuclei calculated on the basis of theoretical (connected by lines) 
and experimental (indicated by markers) GT strength distributions. All rates are plotted relative to those calculated by in the 
shell- model with the KB3G interaction (for which the relative rate thus equals unity). The relative EC rates are shown for two 
combinations of stellar density and temperature: (a) pY e =10 7 g/cm 3 and T = 3xl0 9 K and (b) pY e =10 9 g/cm 3 T = lOxlO 9 
K. 



TABLE II. Average deviations between EC rates calculated based on theoretical GT strength distributions and EC rates based 
on GT strengths extracted from charge-exchange experiments, relative to the experimental values, for two stellar density- 
temperature combinations. Only EC rates on the ground states of 7 nuclei ( 48 Ti, 51 V, 58 Ni, 60 Ni, 62 Ni, 64 Ni and 64 Zn) in the 
pf shell for which high-resolution data were available are considered, or the strength distribution at low excitation energies 
in the daughter nucleus were otherwise well-defined. The left-hand side of the table refers to deviations at pY e =10 7 g/cm 3 
and a temperature of T = 3xl0 9 K (case I) and the right-hand-side of the table refers to deviations at pY e =10 9 g/cm 3 and 
T = 10 x10 s K (case II). 



I: pY e =W 7 g/cm 3 T = 3xl0 9 K II: pY e =10 9 g/cm 3 T = lOxlO 9 K 

GXPFla KB3G QRPA GXPFla KB3G QRPA 

Abc -0.24 -0.34 29. -0.05 0.01 0.66 



A EC \ 0.31 0.47 30. 0.08 0.30 0.72 



experimental GT strength distributions. Average abso- 
lute deviations were less than 50%. At higher densities 
and temperatures (pY e =10 9 g/cm 3 and T = 10xl0 9 ii), 
these deviations became smaller: the average absolute 
deviation between the EC rates based on the shell-model 
calculations using the GXPFla interaction and the data 
was less than what could be attributed to experimental 
uncertainties (10%). For the calculations with the KB3G 
interaction that deviation was slightly higher (30%), but 
when the sign of the deviations were taken into account, 
the average deviation reduced to below the uncertain- 
ties. EC rates calculated on the basis of GT strength 
distributions calculated in QRPA exhibited much larger 
deviations, in particular at low stellar densities. Usage of 
EC rates based on these theoretical calculations is better 
avoided in simulations of late stellar evolution that are 
very sensitive to these rates. 

The results from this work can serve as a benchmark 
for future theoretical studies that aim to provide esti- 



mates for EC rates of relevance for stellar evolution. In 
addition, they give a first estimate for the uncertainties 
in EC rates that are used in astrophysical simulations 
which can serve as input for sensitivity studies in stellar- 
evolution modeling. 



Except for the lightest nuclei (near 45 Sc), sufficient 
data are available to test theoretical predictions for GT 
strength distributions in stable p/-shell nuclei. However, 
very little is known about the quality of the theoretical 
predictions and the derived EC rates for p/-shell nuclei 
further away from the valley of stability, and for nuclei 
beyond the p/-shell, since data are scarce and the avail- 
able theoretical tools more limited. Future experimental 
programs and advances in theoretical techniques will al- 
low for extension of comparative studies such as the one 
presented here to these heavier and more exotic nuclei. 
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